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THE  LINDHARO,  NIELSEN  AND  SCHARFF  METHOD  OF 
OBTAINING  APPROXIMATE  CROSS  SECTIONS 


INTRODUCTION 

Classical  scattering  theory  gives  an  exact  expression  for  generating  a 
differential  cross  section  (6o50),  or,  equivalently,  the  relation  between 
scattering  angle  and  impact  parameter,  from  a  known  interparticle 
potential.  Conversely,  if  one  has  an  exact  cross  section  there  is  an  exact 
inversion  process  that  yields  the  potential  that  is  the  counterpart  of  that 
cross  section. 

Linohard,  Nielsen  and  Scharff  (LNS)  developed  an  approximate  method  to 
generate  a  cross  section  from  a  potential  and  used  this  method  {Li68)  to 
generate  cross  sections  from  both  power  law  potentials  ano  the  Thomas-Fermi 
potential.  Their  Thomas-Fermi  cross  section  is  widely  useo  in  atomic 
scattering  and  radiation  damage  calculations. 

Given  the  approximations  that  make  up  the  LNS  method,  there  is  an 
inversion  procedure  that  generates  the  original  potential  from  the 
approximate  cross  section. 

In  the  following  sections  we  describe  the  three  approximations  that  make 
up  the  LNS  method  and  the  corresponding  inversion  process.  Vie  furnish 
several  examples,  including  one  that  demonstrates  the  utility  of  starting 
with  a  cross  section  and  using  the  simplifiea  inversion  process  to  generate 
the  approximate  potential  corresponding  to  that  cross  section. 

Manuscript  approved  May  21, 1986. 
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EXACT  SCATTERING  RELATIONS 


Consider  a  particle  of  energy  E,  Charge  Z^e  and  mass  interacting 
with  a  stationary  particle  of  charge  Z2e  and  mass  M2  witn  the 
interaction  represented  by  the  potential  V(r).  If  we  write  the  potential  in 
the  form 

V(r)  =  A  C(x)  (2-1 

X  »  r/a, 

where  a  is  some  characteristic  length,  then  the  classical  relation  between 
the  angle  through  which  the  impinging  particle  scatters  and  its  impact 
parameter  is  given  by  (6o50) 

e  =  n  -  2p  f  dr  r"^  [  1  -  _  p2/  ]-l/2  ^  (2-2 

where  p  is  the  impact  parameter,  r^  is  the  distance  of  closest  approach 
(the  largest  root  of  the  radical  in  Eq.  (2-2)),  and  e  is  the  dimensionless 
energy 

€  =  E^  /  A  =  E  [M2/(Mi+M2)]  /  A  =  E  /  E^  ,  (2-3 

where  E  is  the  center  of  mass  energy  corresponding  to  lab  energy  E  and 
the  Lindhard  energy  unit  is 

E^  =  A  (M^+M2)/M2  .  (2-4 

2  1 

(For  the  Coulomb  interaction  l\sZ^z^e  la.  and  C(x)»x“  . ) 

The  relation  between  the  energy  T  transferred  in  a  collision  and  the 
scattering  angle  is 

T  =  T^  sin^  ©/2  ,  (2-5 


where 


^  ”2  . 
'm  *  T  ^ 

(Mj  + 


is  the  maximum  kinetically  allowed  energy  transfer. 
The  differential  cross  section  is  given  by 


da  *  2  w  p  dp  ■  -2*  a^  sin  e  de  a  -oj  dT  . 


LINDHARD-NIELSEN-SCHARFF  SCATTERING  RELATIONS 

From  a  desire  to  simplify  their  atomic  scattering  and  raoiation  damage 
calculations,  Lindhard,  Nielsen  and  Scharff  (LNS)  (Li68)  made  a  series  of 
approximations  co  the  relations  between  the  energy,  impact  parameter, 
scattering  angle  and  cross  section.  They  started  by  invoking  the  momentum 
approximation  'm',  which  gives  the  scattering  angle  as  a  function  of  energy 
and  impact  parameter  to  be  (Le63) 


•m  '  -P'' 


(We  will  use  the  notation  of  identifying  letters  in  single  quotes  to 
indicate  which  approximations  are  being  used.) 

Instead  of  the  energy  transfer  T  they  use  the  reduced  variable 

r?-  =  T/T^  »  sin^  e/2  .  (3-2 

They  then  define  what  they  call  the  wide  angle  approximation  'w',  which 
consists  of  the  using  the  relation 

4'^  ‘3-3 

to  replace  e,  wherever  it  appears,  by  2?i/e.  Upon  using  this  expression  in 
Eq.  (3-1)  we  obtain 


-.v-  •  V- 


,oo 

"mw -- I  1.  ^C(x)  (3-4 

y  p/a 

Equation  (3-4)  provides  a  functional  relationship  between  the  impact 
parameter  and  the  reduced  energy  transfer  that  can  be  used  to  obtain  the 
cross  section.  The  advantage  of  the  LNS  method  is  that  the  differential 
cross  section  depends  on  only  one  variable  (either  t  or  p/a);  ordinarily 
there  would  be  an  explicit  dependence  on  both  p  and  E.  LNS  now  write 

doj^  =  2it  p  dp  =  2»  p(n)  dn  »  -  ir  a^  n  ^  f(n)  dn  .  (3-5 

which  defines  tne  function  f(n)  as 

fi  \  o  2  —2  /  \  dp(n)  2  —2  .  dn(p)  “i”!  /*>  c 

f(n)  =  -  2  n  a  p(n)  ^  =  -  n  a  L  J  •  (3-6 

Me  will  call  f(n)  the  kernel  of  the  cross  section. 

The  last  approximation  is  what  we  might  call  the  kinematical  correction 
'k',  which  insures  that  T»yE  as  p»0.  This  is  not  the  case  with  expression 
(3-4),  which  yields  as  arbitrarily  large  value  of  n  as  p-^O.  We  cure  this 
problem  when  evaluating  Eq.  (3-4)  by  replacing  p  by 

p  -  [p^  ♦  PoJ^'^  (3-7 

where  p^  has  the  value  such  that  n*€  when  p»0.  Thus  the  expression  for  n 
becomes 

WfP)  -  ■  <3-8 

The  approximation  'k'  need  only  be  invoked  when  explicit  use  is  being  maoe 
of  the  impact  parameter.  When  working  with  the  differential  cross  section 
as  a  function  of  n.  the  impact  parameter  does  not  appear. 
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TOTAL  CROSS  SECTION  AND  STOPPING  POWER 


The  total  cross  section  corresponding  to  expression  (2-7)  would  be 


but  the  total  cross  section  is  infinite  for  the  infinite  range  potentials 
that  we  consider.  One  can  avoid  this  problem  by  imposing  a  maximum  allowed 
impact  parameter  or  a  cutoff  on  the  energy  transfer  allowed  at  low 
energies. 

The  stopping  power,  which  is  the  energy  loss  per  unit  path  length,  is 
given  by 


roo 

S(E)  =  2  ir  /  dp  p  T(p) 


f(n)  . 


(4-2 


The  stopping  power  is  finite,  even  for  infinite  range  potentials. 

We  can  use  Eq.  (3-5)  to  obtain  relations  between  the  impact  parameter 
and  the  reduced  energy  transfer;  within  the  spirit  of  the  LNS  approximations 
these  relations  are  the  inverses  of  (3-4)  and  (3-8).  For  the  'mw'  case  we 
have 


and  for  the  'mwk*  case 


(4-3 


q 


2 


(4-4 


The  latter  result  is  the  more  meaningful  physically,  but  we  will  have  need 
of  Eq.  (4-3)  when  we  examine  the  inverse  relations  that  generate  potentials 
from  approximate  cross  sections. 


EXACT  INVERSION 


Just  as  we  can  calculate  the  scattering  angle  from  the  potential  using 
Eq.  (2-2),  there  is  an  inverse  relation  that  yields  the  interaction  when  the 
scattering  angle  is  furnished  as  a  function  of  the  impact  parameter.  Firsov 
(Fi53,  To72)  established  that  if  we  define  a.  as 

J  [1  -  V(r)/E^j  (5-1 

where  E^  is  the  center  of  mass  energy  (Eq.  (2-3)),  then  the  interaction 
separation  x  is  related  to  o(p)  by 

^00 

x(a))  =  0)  expL  -^  /  dp - i;2  J 

6y  using  Eqs.  (5-1,2)  one  can  obtain  V(r).  If  one  uses  the  exact  inversion 
procedure  on  an  approximate  cross  section,  the  resulting  interaction  will  be 
energy  oepenoent.  Robinson  (Ro69)  shows  this  effect  using  the  LNS  cross 
section. 


IMPACT  expansion  INVERSION 


The  equivalent  to  the  Firsov  procedure  when  one  is  using  the  momentum 
approximation  'm'  to  relate  the  scattering  angle  to  the  energy  and  impact 
parameter  is  the  expression  (Sm66,  To72) 


2  r  ®m^P^ 


00 


(6-1 


In  reducea  form,  this  becomes 


C(x; 


2 

—  e 
ir 


2 

—  e 

w 


/. 

L 


CO 

dp 


dp 


(6-2 


■■  V 


.  •  .  •  s.'  *  I 


and  upon  invoking  the  wioe  angle  approximation  we  nave 


C(x) 


4 

IT 


(6-3 


One  can  easily  show  that  Eq.  (3-4)  ana  Eq.  (6-3)  form  an  inversion  pair  (by 
using  Eq.  (3-4)  to  substitute  for  in  Eq.  (6-3),  interchanging  the 
oroer  of  integration,  and  evaluating  the  resulting  integral),  as  do  Eq. 
(3-1)  and  Eq.  (6-2).  Finally,  if  one  has  the  'mwk'  version  of  n,  one  uses 


C(x) 


^mwk 


(6-4 


One  should  be  consistent  in  using  these  approximations  when  inverting  to 
obtain  a  potential . 


AN  EXAMPLE:  THE  INVERSE  SQUARE  POTENTIAL 

Let  us  assume  an  inverse  square  scattering  potential,  given  by 

V(r)  =  A2  x"^  .  (7-1 

Eq.  (2-2)  yields  the  exact  scattering  angle  and  reduced  energy  transfer 

e  =  It  [1  -  q(q2  +  l/e)"^^^j  (7-2 

n  =  e  cos  irq/(2(q^+l/€))  . 

In  the  momentum  approximation  the  energy  transfer  becomes 

Ti|^  =  e  sin  Tr/(4cq^)  .  (7-3 


With  the  the  use  of  the  wide  angle  ana  kinematic  approximations  the  reduced 
energy  transfer  becomes 


and 


W  =  ^/(l-"4eq2/ir)  ,  (7-5 

respectively. 

We  can  compare  the  various  approximations  for  tne  inverse  square 
potential  Oy  evaluating  the  exact  potential  that  corresponds  to  each  of  the 
approximate  energy  transfer  relations,  that  is,  to  each  approximate  cross 
section.  We  use  the  exact  inversion  in  the  form  of  Eq.  (5-1)  ana 

Jroo  _i 

dq  j  •  (7-6 

u)  (q  -  «  )  ' 

2 

Of  course,  if  we  use  the  exact  expression  (7-2)  we  obtain  C(x)=x  .  For 

2 

the  others,  we  can  make  expansions  in  1/x  ,  thus  obtaining  the  results  in 
the  form 

msO 

The  first  few  coefficients  b^^  for  the  three  approximations  are  given  in 
Table  1. 


Table  1 


Coefficient  b^,  of  e"'^  x-2m-2  in  Eq.  (7-7) 


Inaex 

Exact 

Value 

'm' 

Approximation 

'mw' 

'mwk ' 

0 

1 

1 

1 

1 

1 

0 

.500U0 

.50000 

-.02360 

2 

0 

.66667 

.72150 

.00329 

3 

0 

1.12500 

1.28949 

-.00071 

4 

0 

2.13333 

2.58358 

.00019 
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We  see  that  the  comoined  'mwk'  approximation  gives  a  cross  section  whose 
corresponding  potential  more  closely  agrees  with  the  exact  cross  section 
than  do  the  other  approximations.  In  this  regara  we  should  mention  that  for 
the  Coulomb  potential,  the  'mwk'  approximation  is  identical  to  the  exact 
cross  section. 


AN  EXAMPLE:  THE  POWER  LAW  KERNEL 

There  are  cases  in  which  we  start  with  a  cross  section  directly,  rather 
than  calculating  one  from  a  potential.  Consider  the  kernel 


=  B  n 


1-2 /n 


(8-1 


By  using  the  relation  (4-3)  we  have 


■r 


ds  s~^  B 


n  c  _-2/n 

7 


B  n 


(8-2 


or 


j-n  B  Tn/2 

n  =  L - 

2 


(8-3 


as  the  relation  between  the  impact  parameter  and  reduced  energy  transfer. 
Upon  using  the  inversion  formula  (6-3)  we  fine 


r/xl  4<nB^n/2  • 


n-1 

dB  sin"  ^6 


(8-4 


The  various  constants,  including  the  integral,  don't  interest  us,  but  we  see 
that  our  kernel  is  equivalent  to  a  power  law  x”*^  potential.  We  note  that 
for  the  Coulomo  potential,  for  which  C2^(x)=x~^,  B  will  have  the  value 
B=l/2  and  f^ (n)=l/ (2n) . 

The  point  of  this  exercise  is  to  illustrate  that  we  can  start  with  some 
simple  form  for  the  cross  section  and  use  the  inversion  relations  to 
determine  the  equivalent  potential,  within  the  LNS  approximations. 


AN  EXAMPLE:  A  SIMPLE  KERNEL  FOR  A  FINITE  RANGE  POTENTIAL 


Assume  the  kernel 

fin)  =  (1/2)  n^/(S+n)^  .  (9-1 

We  see  that  for  large  values  of  n,  corresponding  to  nigh  momentum  transfers, 
this  f(n)  tends  toward  the  Coulomb  kernel.  By  using  the  relation  (4-3)  we 
have 

=  (1/4)  (e+n)"^  (9-2 

or 

n  =  l/2q  -  B  (&_3 

as  the  relation  between  the  impact  parameter  and  reduced  energy  transfer. 

We  see  that  zero  energy  transfer  (n=0)  corresponds  to  a  maximum  value  of  the 
impact  parameter  of 

=  1/(26)  .  (9-4 

Upon  using  the  inversion  formula  (6-3)  we  find 

C(x)  =  x“^  [2/Tr]  [cos~^  y  -  y  ln((l+(l-y^)^^^)/y)]  (9-5 

where 


y=x/x^  .  (9-6 

Note  that  upon  setting  b  to  zero  (y=0)  one  recovers  the  reduced  Coulomb 
potential  C(x)=l/x.  The  quantities  in  square  brackets  in  Eq.  (9-5) 
constitute  the  finite  range  screening  factor  for  the  Coulonm  interaction. 

In  Fig.  1  we  illustrate  the  effect  of  this  screening  factor.  Rather 
than  set  b  to  different  values  directly,  we  will  choose  various  values  of 
X|^*1/(2b).  For  atomic  scattering  problems  a  typical  screening  radius 


(Eq.  (2-1))  is  a»0.01  nm,  so  that  Xjj|®100  corresponds  to  a  cutoff  radius  of 
1.0  nm.  Figure  1  shows  the  screening  factor  for  three  potentials  derived 
from  the  cross  section  (9-1),  using  values  of  5,  30,  and  100. 

If  one  had  a  simple,  finite  range,  pnenomenological  potential  with 
Coulomb  behavior  at  small  separations,  one  could  fit  that  potential  by 
adjusting  a  (or  x^^^),  obtaining  in  the  process  a  simple  cross  section 
corresponding  to  that  potential.  This  same  procedure  can  be  used  to  fit  a 
wioe  variety  of  phenomenological  potentials  with  forms  whose  corresponding 
cross  sections  are  simple  (Mu85). 


NUMERICAL  EVALUATION 

When  using  numerical  methods  to  evaluate  the  inversion  integrals,  it  is 
advantageous  to  modify  formula  (6-3)  by  a  change  of  variable  (from  q  to  n) 
and  an  integration  by  parts.  This  yields 


i  h  dn  (q2  -  x2)l/2  A  ,  ^3 

"Jo  ^  r 


(10-1 


where  q  and  q'  are  functions  of  n.  The  variable  is  implicitly  given  by 


dn  n  ^  f(n)  . 


(10-2 


and  the  function  q(n)  is 


os  s~^  f(s)  ; 


(10-3 


q'  is  simply 


q'  .  -n"-"  f(n)  /  (2q) 


(10-4 


so  that 


C(x) 


(10-5 


.i  h  dn  (q2-x2)l/2  ll.  114] 

Jo  ^  2nq 

The  inversion  proceoure  consists  of: 

a)  Choosing  some  value  for  and  using  Eq.  (10-2)  to  determine 
the  corresponding  value  of  x. 

b)  Choosing  a  set  of  integration  points  for  the  numerical 
integration  of  Eq.  (10-5)  and  using  Eq.  (10-3)  to  determine  the 
corresponding  value  of  q. 

Often,  if  the  kernel  f(n)  is  some  simple  form,  then  the  integral  (10-3) 
can  be  done  exactly.  If  Eq.  (10-3)  must  be  evaluated  numerically,  then 
one  can  use  the  form 


DISCUSSION 

To  extend  the  discussion  of  the  last  two  sections,  many  radiation  damage 
calculations,  such  as  Boltzmann  transport  theory  calculations,  use  a  cross 
section  directly  without  any  reference  to  the  potential.  For  such  numerical 
work  a  simple  cross  section  can  save  computer  time. 

Our  earlier  treatment  of  the  Moliere  potential  serves  as  an  example. 

The  Moliere  is  an  approximation  to  the  Thomas-Fermi  interaction  that  has  an 
exponential  tail  instead  of  the  Thomas-Fermi ' s  x~^  dependence  at  large 
separations  (Mo49,  To72a,  Ro69). 

In  order  to  generate  a  useful  cross  section  we  (Mu78,  Mu80): 

a)  Used  the  LNS  procedure  to  generate  the  kernel  f|vi(n) 
corresponding  to  the  Moliere  potential. 

b)  Created  an  approximate  fit  to  f|^(n)  that  was  expeditious  for 
machine  calculations. 

c)  Used  the  approximate  inversion  procedure  to  examine  the  aeviation 
(within  the  LNS  approximations)  of  the  potential  corresponding  to 
our  fit  and  the  Moliere  potential. 


This  method  allowed  us  to  find  a  simple  form  for  the  cross  section  that 
reproduced  the  Moliere  potential  to  within  six  per  cent. 

The  simplicity  of  the  LNS  inversion  pairs  generally  makes  it  easy 
to  create  a  usable  fit  to  the  cross  section  corresponding  to  a  given 
atom-atom  interaction. 
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